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Abstract 



This paper analyzes the Lamb vector divergence, also called the hydrody- 
namic charge density, and its implications to the Navier-Stokes system. It 
is shown that the pressure field can be always chosen in a way that ensures 
regularity of the Navier-Stokes system. The abstract pressure field that en- 
sures regularity is defined through two partial differential equations, being of 
the elliptic kind. The pressure field defined such a way can be interpreted as 
a control potential field that keeps the system regular. The controlling pres- 
sure field depends only on the velocity field of the fluid and its derivatives, 
so that the result is applicable in any general setting where the initial data 
is divergence free, smooth and square-integrable. 



1 Introduction 



The Navier-Stokes equations are supposed to model any kind of fluid flow, 
including turbulent ones. The essential problem with the equations is the 
nonlinearity of the convective term. The pressure gradient represents trans- 
port and the diffusive term represents loss of energy in the system. In other 
words, to understand the equations, one has to understand the nonlinear 
convection term. The Navier-Stokes equations can be written in vector form 
as follows 

9u _ . _ _ 

— — h u • vu = uAu — \p (1) 
at 

where u(x,y, z,t) G M. 3 is the velocity field of the fluid and p(x,y,z,t) G M 
is the scalar pressure field. The diffusive term depends on a constant v > 0. 
The convective term is u ■ Vu. This term can be decomposed according to 
the vector calculus identities 

u ■ Vu = u x u + V(^m- u) (2) 

Subsituting this into the Navier-Stokes equations, one has 

du , 1 _ _ , . 

— = — u x u — V(jp + -u ■ u) + pAu (3) 

where u = V x u is the vorticity of the fluid flow. This formulation of the 
Navier-Stokes equations is known as the Lamb formulation [lj. 

In addition we assume that the fluid flow is incompressible, that is 

V ■ u = (4) 

The system is assumed to have some smooth and square-integrable initial 
data u° = u(x, y, z, 0) over the whole space M 3 . 



2 The regularity of Navier-Stokes equations 

When we say that the solutions are regular, we mean that two conditions are 
fulfilled, see 0. 

u,p G C°°(M 3 x [0,oo)) (5) 

and 

| u | 2 dx < C (6) 
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for all t > 0. It is a well known fact that the solutions are regular if the 
enstrophy of the system stays bounded, see [3J. This means that we will 
require 

u | 2 dx < c (7) 

for all t > 0. This implies that it is sufficient that the vorticity does not blow 
up. 



3 The Lamb vector and its divergence 

The Lamb vector L is the nonlinear term in the Lamb formulation of the 
Navier- Stokes equations 

L = u x u (8) 

By applying the divergence operator to the Navier-Stokes equations, one has 
for the divergence of the Lamb vector 

V • (uj x u) = -A(p+ hi- u) (9) 

Where the Laplacian is the normal, or scalar Laplacian. Marmanis and Shri- 
dar have called the divergence of the Lamb vector as hydrodynamic charge 
density, see [5], [B] in analogy to the electromagnetic theory. 



3.1 The divergence theorem and Lamb vector diver- 
gence 

The Lamb vector divergence can be represented as a sum of two parts 

V ■ L = u ■ V x u — u ■ u (10) 

,see an excellent paper on Lamb vector divergence by Hamman et al |7J. The 
first term, u ■ V x oj is called the flexion product whereas the second term 
is related to enstrophy. From the vector calculus, we have the divergence 
theorem: 

tf](V-F)dV = §F.S (11) 

v s 
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where V is a compact subset of M 3 and F is some continuously differentiable 
vector field in V. The vector S is the normal vector to a smooth surface 
embedding V. Moreover, from the vector calculus identities we have 

JJJ (V • L)dV = JJ[(w • V x Cj - u ■ ou)dV = (jj(u; x u) ■ S (12) 

V V s 

Now let us substitute the explicit Lamb divergence on the left side of the 
equation 

JJ[(-A(p + ■ u))dV = JJJ (w • V x u - u ■ u)dV (13) 

V V 

We note immediately that we can obtain the local enstrophy from the equa- 
tion 



JJJ uj ■ udV = JJJ(w • V x u)dV + JJJ(A(P + if ■ u))dV (14) 



V V V 

Putting the right side under the same integral 



JJJw-dW = JJ[(£- V xu + A(p+^u-u))dV (15) 



V V 



4 Pressure as stabilizing force ensuring reg- 
ularity 

Now, in order to prevent a blow-up, a sufficient condition is to have the 
enstrophy bounded in IR 3 , so we want the left side of the equation to be 
bounded from above. We will proceed in two steps. First, let us suppose 
now that the pressure field p is composed of two parts as follows 

p(x, y, z, t) = q(x, y, z, t) + r(x, y, z, t) (16) 

Substituting this into the integral equation yields 

JJJ u ■ udV = JJJ(m -V xu + A{q + r + ^u- u))dV (17) 
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Taking the Laplacian term by term 

JJJ u ■ cudV = JjJ(« -V xu + Aq + Ar + A(^u ■ u))dV (18) 

R 3 R 3 

We want to ensure that the first component of the pressure field partly bal- 
ances the stress from the convection terms by setting 

Ag = -u- V x u - A(]-u ■ u) (19) 

This pressure field q is well defined through the Poisson-like equation. What 
is left is 

JJ[ u ■ udV = J j J ArdV (20) 

R 3 R 3 

Differentiating with respect to time we have, 

jM a - adV= jM AriV < 21 » 

R 3 R 3 

This is the control equation for the evolution of enstrophy. We now require 
that this satisfies 

|J0W<O (22) 

R 3 

There is at least one expression that makes this possible and that depends 
only on the velocity field, namely the time evolution of the total kinetic 
energy K(t) 



dK(t) _ d 1 
dt ~ dt2 

So let us choose 



u\ dV <0 (23) 



1 , -.,2 



Ar = - | u | 2 (24) 



Now we have 

d err _ _ dl 



dt 



{{{u udV = ~ [ \u\ 2 dV<0 (25) 
JJJ dt 2 



That is, the evolution of enstrophy is equal to the evolution of total kinetic 
energy. To take stock, we have a pressure field p such that 

p(x, y, z, t) = q(x, y, z, t) + r(x, y, z, t) (26) 
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such that the pressure components satisfy 



Ag = -«-Vxw-A(^5-M) (27) 



and 



At- = - | u T (28) 



Putting them together 



Ap = A(q + r) = -u ■ V xu- A(-u ■ u) + - | u | 2 (29) 

This Poisson equation for the pressure ensures that the time evolution of 
enstrophy satisfies 

which in turn implies regularity for the Navier-Stokes system with smooth, 
divergence-free and square-integrable initial data. 



5 Consequences on the energy dissipation rate 

It is a well known fact, see [3] that the following equality holds 

f = -„ E (3!) 

where K(t) = | J R3 \ u | 2 dV is the total kinetic energy of the system and 
E{t) = J R 3 | u | 2 dV is the total enstrophy of the system. Let us differentiate 
the energy evolution equation with respect to time 

d 2 K(t) dE(t) 

and by noting that dK Jf' = dE Jjr we will have the result 

d 2 K{t) dK(t) 

which has the general solution 

K{t) = e rt (34) 
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where r = or r = —v. The case where r = refers to the specific case 
where the initial data has zero vorticity. To exclude that trivial case, we 
assume r = — v and thus have 



so that the total kinetic energy of the system and thus vorticity decay expo- 
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